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ABSTRACT 

Context. Galaxy clusters can be used as cosmological probes, but to this end, they need to be thoroughly understood. Combining all 
cluster observables in a consistent way will help us to understand their global properties and their internal structure. 

Aims. We provide proof of the concept that the projected gravitational potential of galaxy clusters can directly be reconstructed from 
X-ray observations. We also show that this joint analysis can be used to locally test the validity of the equilibrium assumptions in 
galaxy clusters. 

Methods. We used a newly developed reconstruction method, based on Richardson-Lucy deprojection, that allows reconstructing 
projected gravitational potentials of galaxy clusters directly from X-ray observations. We applied this algorithm to the well-studied 
cluster Abell 1689 and compared the gravitational potential reconstructed from X-ray observables to the potential obtained from 
gravitational lensing measurements. We also compared the X-ray deprojected profiles obtained by the Richardson-Lucy deprojection 
algorithm with the findings from the more conventional onion-peeling technique. 

Results. Assuming spherical symmetry and hydrostatic equilibrium, the potentials recovered from gravitational lensing and from X- 
ray emission agree very well beyond 500 kpc. Owing to the fact that the Richardson-Lucy deprojection algorithm allows deprojecting 
each line of sight independently, this result may indicate that non-gravitational effects and/or asphericity are strong in the central 
regions of the clusters. 

Conclusions. We demonstrate the robustness of the potential reconstruction method based on the Richardson-Lucy deprojection 
algorithm and show that gravitational lensing and X-ray emission lead to consistent gravitational potentials. Our results illustrate the 
power of combining galaxy-cluster observables in a single, non-parametric, joint reconstruction of consistent cluster potentials that 
can be used to locally constrain the physical state of the gas. 

Key words, galaxies: cluster: general - X-rays: galaxies: clusters - gravitational lensing: strong - gravitational lensing: weak 


1. Introduction 

The cold dark matter (CDM) model favor ed by substantial theo- 
retical and observational evidence (e.g., iKomatsu et al.luOllh 
predicts the formation of small structures at early times and 
their later merging into larger systems. In the CDM framework, 
galaxy clusters are expected to be the latest gravitationally bound 
systems to form. 

Galaxy clusters are mainly composed of three components: 
galaxies, hot gas, and dark matter. The dark matter is by defini¬ 
tion impossible to observe directly, but it is expected to follow 
a universal density pr ofile approximated b y the Navarro-Frenk- 
White profile (NEW. iNavarro et'aDll997l) . However, observa¬ 
tions of the intracluster medium (ICM) and of cluster galaxies 
can constrain the dark-matter distribution itself. As the ICM fills 
the deep gravitational potentials, its distribution is governed by 
the depth and sha pe of the dark-matter potential well (see, e.g., 
lEttori et ^120131 for a review), hence the observables provided 
by the thermal plasma carry information on the gravitational po¬ 
tential. 

The physical properties of the ICM can be directly studied in 
X-rays through thermal bremsstrahlung and line emission and at 
millimetrer wavelengths through the Sunyaev-Zel’dovich (SZ) 
effect. On the other hand, measurements of gravitational lensing 
in clusters constrain the gravitational tidal field and thus the cur¬ 


vature of the projected gravitational potential of these clusters. 
Our study is based on the idea that a simultaneous combination 
of as many galaxy-cluster observables as possible may lead to a 
consistent model for the projected cluster potential because un¬ 
der equilibrium assumptions, all these observables can be com¬ 
bined on the grounds of the gravitational potential. 

In addition to improving the reconstruction of the projected 
cluster potential beyond gravitational lensing, combining as 
many cluster observables as possible will also allow us to cover 
a wider range of angular or radial scales. The joint analysis 
of weak and strong gravitational lensing allows us to map the 
dark-matter distribution within the entire cluster. Furthermore, 
because the emissivity depends on the square of the gas den¬ 
sity, X-ray flux and temperature profiles ar e most reliable in 
and near the central region of the cluster (lEttori & Molendil 
1201 ih . while the thermal SZ effect has a shallower dependence 
on the local gas density because the thermal SZ signal is pro- 
portional to the gas pressure integrated over the line of sight 
(iPlanck Collaboration et alJ 1201 3t Eckert et"^ l2013h . Hence, 
owing to the different dependences of different observable sig¬ 
nals on angular scales, weak-lensing and thermal-SZ measure¬ 
ments extend to relatively large distances from the cluster cen¬ 
ter, while strong-lensing and X-ray observations are more sen¬ 
sitive to regions near the cluster center. In combination, more 
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reliable reconstructions of cluster potentials can be expected on 
all scales. 

Another advantage of this joint analysis is that the assump¬ 
tions made on the physical state of the gas can be verihed. The 
lensing potential can be recovered without any equilibrium as¬ 
sumptions, while the gravitational potential reconstructed from 
X-ray and thermal-SZ observations is recovered assuming hy¬ 
drostatic equilibrium. Therefore, the projected gravitational po¬ 
tential recovered from gravitational lensing observations can be 
compared with the gravitational potential recovered from the hot 
ICM observations, and the result can be used to determine which 
clusters satisfy the assumptions and which do not. This will help 
us to distinguish the different processes occurring in clusters. 

We here focus on reconstructing the projected potential 
of galaxy clusters using X-ray observations. However, other 
observables may also be used, such as galaxy kinematic s 
(ISarli et al.ll2013l) or the thermal SZ effect (IMaier et al.ll2013h . 
Our hnal goal is to set up one common likelihood function for 
the gravitational potential, joining information from all available 
observables. As a proof of concept, in this study the projected 
potential of the well-studied cluster Abell 1689 is reconstructed 
from X-ray observations. 


The cluster Abell 1689 is located at redshift 0.183 in the 
Virgo constellation. As a well-known strong-lensing cluster, it 
has been the subject of many studies. The multiwavelength 
observations performed on this object extend from the mid- 
infrared, for instan ce, using ISOCAM onboard the ISO satellite 
( Fadda_et^lj2000h, to the X-ray band, for ins t ance, u sing Suzaku 
( Kuwharada_e^ ] l2010h . ROSAT (lAllenI T 998h. Chandri^ 
( Xue &WtJ |2QQ2l)! and XM M-Newton (lAndersson & Madeisa 
l2004t ISnowden et al.ll2008h . In visible light, the optical signal 
has been observed using the Subaru/Suprime-Cam te lescope and 
the Hubble Space telescope (iBroadhurst et al.ir2005h . This list is 
not comprehensive. 

The main result of these observations is that Abell 1689 
is a dynamically active cluster: kinematics and X-ray stud¬ 
ies i ndicate a merger aligned wi th the line of sight (see, 
e.g., lAndersson & Madeiskil |2004 and references therein). 
Combinations of different observations in joint analyses have 
been performed by many authors and have shown that the hy¬ 
drostatic mass is lower than the mass estimated from gravita- 
tional lensing measurements in the central regions (see, e.g., 
Peng et ^ 2009t Morandi et alJ 20 ill Miralda-Escude & Babui 


IHIl Lokas et al Jl2006HSereno et al.l2013h. 

iMiralda-Escude & Babull (Il995h combined strong-lensing 
and X-ray surface brightness measurements to test the physi¬ 
cal state of the gas near the cluster center. The authors showed 
that the observed temperature in the central regions is lower than 
that expected from hydrostatic equilibrium. As a result, a non- 
negligible nonthermal pressure component should be present 
and act against gravity. According to these authors, this nonther¬ 
mal pressure support could be ass ociated with the m erger pro¬ 
cess. In a complementary analvsis. lLokas et ^ (l2006h measured 
the velocity distribution of the galaxies in the held and identihed 
several matter clumps aligned with the line of sight through the 
cluster, which do not interact with the cluster dynamically and 
might affect the lensing mass estimates without modifying the 
X-ray mass estimate. This provides an alternative explanation 
for the mass discrepancy observed be tween the two me thods. 
In a different analysis of Abell 1689, iPeng et aTI (l2009h com- 


; also the Chandra catalog of ICavagnolo et all ll2009h . 


bined weak-lensing with X-ray data from the Chandra satellite. 
Considering the central part of the cluster out toR ~ 1000 kpc/h, 
they found that at radii larger than R ~ 200 kpc/h, the X-ray and 
lensing masses are consistent, whereas they found a discrepancy 
between the masses at smaller radii. These authors explained this 
discrepancy by a projection effec t. This conclusion agrees with 
the results of ISereno et al.l (l2013h . In this study, the authors con¬ 
strained the shape and orientation of the cluster in a joint analysis 
of the ICM distribution in the cluster. Using X-ray and SZ mea¬ 
surements, the cluster appears to have a triaxial shape elongated 
along the line of sight. According to these authors, this may ex¬ 
plain the mass discrepancy. 

Finally, in a joint analysis of X-ray s urface brightne s s with 
strong- and weak-lensing measurements, iMorandi et al.l (1201 ih 
inferred the fraction of nonthermal pressure required for the 
cluster to be in hydrostatic equilibrium, taking into account the 
triaxial shape of the cluster. The authors found that 20 % of the 
total pressure should be nonthermal, most probably contributed 
by turbulent gas motion driven by hierarchical merger events. 

To summarize, a mass discrepancy has been observed in 
many studies, and different explanations have been suggested. 
To understand the origin of this discrepancy, we here propose 
the following method: 1) reconstruct the three-dimensional gas 
distribution using a deprojection method that operates indepen¬ 
dently on each l ine of sight: the R ichardson-Lucy deprojection 
algorithm fR-L. ILuc^I 1974 [19941) : and 2) substitute the mass 
with a local quantity, that is, the gravitational potential. 

This procedure, applied here to a joint analysis of X-ray and 
lensing observations, is expected to allow us to directly compare 
the dark matter and gas distributions and to test the validity of 
the equilibrium assumptions at each projected radius. Using the 
gravitational potential should render the comparison with gravi¬ 
tational lensing measurement more straightforward because the 
gravitational potential is directly observable from the lensing 
measurements. 

To recover the cluster potential from X-ray measurements, 
we hrst deproject the X-ray data, then convert the X-ray three- 
dimensional prohle into the three-dimensional gravitational po¬ 
tential using equilibrium assumptions, and project it again to 
compare it with the two-dimensional gravitational potential ob¬ 
tained by gravitational lensing. Applying this method with dif¬ 
ferent geometries (spherical and triaxial) is expected to remove 
the degeneracy between the different possible explanations for 
the mass discrepancy. We test this algorithm here for the first 
time on a real cluster assuming spherical symmetry. The gener¬ 
alization to a spheroidal shape is ongoing. 

This paper is structured as follows: In Sect. |2] we list the 
assumptions of the method. In Sect. [3] we apply the R-L de¬ 
projection algorithm on the cluster Abell 1689 and compare 
this with the onion-peeling deprojec tion method (O-P, see, e.g., 
iFabian et '^119811: iKriss et al.ll 19831 for a detailed overview of 
this method), which is the standard deprojection method in X- 
ray astronomy. In Sect. 13.11 O-P and R-F deprojection methods 
are directly compared for different cluster observables, and in 
Sect. 13.21 the advantage of the R-F method with respect to the 
O-P method is discussed. In Sect. 13.31 the reconstructed pro¬ 
jected gravitational potentia l obtained using the method devel¬ 
oped in iKonrad et al.l (1201 3h is compared to the lensing potential 
of the cluster as recovered by J. Merten (private communication). 
Finally, the results are discussed in Sect. |4] and we conclude in 
Sect.|5] 
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2. Assumptions of the method 

The ge neral reconstruction formalism outlined in iKonrad et alJ 
(1201 3h is based on the following assumptions; 

- The ICM is in hydrostatic equilibrium, 

VP = -pV<D , (1) 

where O denotes the Newtonian gravitational potential, 
while P and p are the gas pressure and density. 

- The gas follows the poly tropic relation 


in the 0.4 - 2 k eV band. For the de tails of the analysis proce¬ 
dure, we refer to lEckert et al.l (1201 ih . 

In this section, we aim at comparing the results of the R-L 
deprojection method with those obtained with the O-P method. 
For this purpose, we extracted the emissivity, the density profile, 
and the temperature profile from the X-ray surface brightness 
using both methods. 


3. 1. Testing the R-L deprojection method 



where the suffix 0 denotes fiducial values for the pressure 
and for the gas density at an arbitrary radius ro and where y 
is the polytropic index. 

- The gas is ideal, 

P = , (3) 

m 

where T is the gas temperature, is Boltzmann’s constant, 
and m is the mean mass of a gas particle. 

- The bremsstrahlung emission is assumed to dominate the X- 
ray observations. This assumption is justified as long as the 
temperature of the plasma does not fall below ~ 2 keV and 
the Fe XXV-Fe XXVI lin e complex at 6 - 7 keV is avoided 
dPeterson & Fabianll2006l) . The bolometric bremsstrahlung 
emissivity, j^, can be written as 

A = CT^^V , (4) 


where C is an amplitude formed by all relevant physical con¬ 
stants. We integrated the emissivity over the energy band 
[0.4 - 2] keV ( 7 [o. 4 ; 2 ]) instead of over the entire electromag¬ 
netic spectrum. Introducing finite boundaries for the energy 
bands modifies the temperature dependence of the emissiv¬ 
ity, but not its dependence on density. Since the density pro¬ 
file is steeper than the temperature profile, and the emissivity 
depends on the square of the density, the emissivity profile is 
dominated by the density profile and uncertainties in the tem¬ 
perature dependence can be considered negligible for now. 
Except for the natural constants and the averaged parame¬ 
ters, all quantities introduced here are assumed to depend on 
the cluster-centric radius. This dependence is not written ex¬ 
plicitly to simplify the notation. The effective adiabatic index 
is assumed to be constant across the cluster. While this has 
been shown not to be true in general, it can be a good ap- 
pro ximation as long as the core of the cluster is avoided (see. 


e.g,. Tozzi & Normanll2001l : ICapelo et al.ll2012L lEckert et alJ 

I2013I and the discussion section). 

The cluster is assumed to be spherically symmetric. 


3. Deprojection of the cluster Abell 1689 


To construct the surface-brightness profile of the cluster, 
we used an archival 13.4 ks observation of Abell 1689 
with ROSATfPSPC. This choice was motivated by the low 
background of this instrument, which allowed us to ex- 
tract high-quality surface-brightness pro files out to large radii 
dVikhlinin et al.lll999l: lEcke rt et al.l 201 Ih. W e reduced the data 
using the ESAS software ( Snowden et al.l n~994ll . The non-X- 
ray background was modeled and subtracted from the observed 
image, and a vignetting correction was applied. The surface- 
brightness profile was then extracted from the corrected image 


The R-L algorithm is an iterative method that tends to depro¬ 
ject the features of the observed projected function on all scales. 
To avoid reproducing the instrumental fluctuatio ns and other 
sources of noise in the deprojected profile, Lucy (lLucvl[l99^ 
suggested to add a regularization term to a likelihood function 
whose scale and amplitude were controlled by two parameters, 
the regularization p arameter a and the smoothing scale L (see 
iKonrad et al.ll2013l for more details). In the following, we adapt 
these parameters to the X-ray data obtained by traditional tech¬ 
niques. 


3.1.1. Three-dimensional emissivity 

The mean value and the error bars of the recovered emissivity 
were obtained using a Monte Carlo (MC) method based on the 
R-L deprojection method as follows; The value of the surface 
brightness profile was randomized within each radial ring, ac¬ 
cording to a Gaussian distribution whose mean was the observed 
value and whose width was given by the error bars. After 10000 
simulations, the deprojected value and its error are given by the 
mean and Icr deviation of the resulting distribution. 

In Fig.[T] we show the results of the X-ray surface-brightness 
deprojection using the O-P and R-L methods. To obtain these 
results, the parameters of the R-L method were set to a = 0.01 
and L - 300 kp c. These parameter values are lower than the 
expectations (see lKonrad et al.l2013l for more details) and imply 
that the regulation term introduced in the R-L method is not very 
effective in the case of Abell 1689. 

More detailed investigations with simulated clusters have 
shown that the results of the reconstruction depend very little 
on the values chosen for the regularization parameters. Testing 
wide ranges for both a and L, typically covering an order of 
magnitude, changes the recovered potentials routinely by much 
less than the statistical uncertainties. 

We emphasize that the regularization term is controlled by 
the difference of the potential reconstructions between subse¬ 
quent iterations steps, which vanishes progressively as the itera¬ 
tion proceeds. The regularization is thus gradually switched off 
toward the end of the iteration. We fixed the number of iterations 
to 10, which corresponds to the number of iterations at which the 
algorithm converges, because no significant changes in the emis¬ 
sivity profile were observed when more iterations were used. 

This is endorsed by our result, which shows that the emis¬ 
sivity profiles obtained with the two methods agree very well 
within the error bars. However, the outermost points of the pro¬ 
file are different in the two techniques, whi ch can be attributed to 
edge effects (see, e.g.. lMcLaughlinlfl999h . The large error bars 
in the outermost points are caused by the weakness of the signal 
at large radii. 
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Fig. 1. Emissivity obt ained by applying the R-L (this analysis) 
and the O-P methods (lEckert et al.ll20i:^ to the ROSAT surface- 
brightness profile. The reconstructed emissivity using the R-L 
method has been normalized to the O-P data points within a 
spherical shell delimited by the radii in the range [183; 2700] 
kpc. 
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3.1.2. Three-dimensional density profile 

The tempera ture of the cluster Abell 1689 lies between 2keV 
and lOkeV (iKawaharada et al.]l2010l) . This means that the soft 
energy band [0.4; 2] keV is below the bremsstrahlung cutoff en¬ 
ergy, and the emissivity only weakly depends on the tempera¬ 
ture. In this particular case, the density profile can be expressed 
as a function of the emissivity, 7 [o. 4 ; 2 ], as 


P(r) oc ^Jj[QA■2]ir) ■ (5) 

In Pig. 12 we compare the density profile obtained using Eq. 13 
for the two depro jection methods: the O-P method (in blue, 
lEckert et al.ll2012l) and the R-L method (in red, this analysis). 
Again, the density profile is well reproduced using the R-L pro¬ 
cedure. 

3.1.3. Temperature profile 

Using the polytropic relation (Eq. Gl) and the fact that we con¬ 
sidered only the soft energy band [0.4;2] keV, the temperature 
profile can be expressed as 

T(r) oc p(ry^' oc ;[o.4;2](r-)‘^‘'^^^ ■ (6) 

Pigure [2 shows the temperature profile reconstructed with a 
polytropic index of y = 1.19 + 0.04 and the temperature pro¬ 
files obtained using different techniques. The polytropic index 
■ y - 1.19 + 0.04 results from fittin g the observed pressure profile 
( Planck Collaboration et al.l2013h with the reconstructed density 
profile (obtained in Sect. 3.1.2), using Eq. l2]i . 

Por completeness, we also fitted Eq. 0 to the joint XMM- 
Newton and Suzaku data and obtained y = 1.17 + 0.03 for the 
polytropic index with Icr error, which is consistent with our ear¬ 
lier result. 

The decrease in the temperature profile supports the validity 
of the polytropic assumption (Eq. 121). 
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Fig. 2. Electron-density reconstructed profile obtained from 
the R-L deprojection method (this analysis) compared to the 
electron- density profile rec overed from the O-P deprojection 
method (lEckert et al.l l2012h . Por both deprojection methods, 
Eq. ||5l has been used. The reconstructed density using the R- 
L method has been normalized to the O-P data points within a 
spherical shell delimited by the radii in the range [183; 2515] 
kpc. 
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3.2. Advantages of the R-L method compared to the O-P 
method 

In observations, the only information available is integrated 
along the line of sight. In the framework of this study, we define 
as local a deprojection method that allows independently depro- 
jecting this projected information for each line of sight, meaning 
that it is not necessary to know the value of the deprojected quan¬ 
tity in one line of sight to derive the deprojected value in another 
and that regions of different projected radius are not mixed dur¬ 
ing the deprojection procedure. 

In Sect. l3.ll we saw that the R-L and O-P deprojection meth¬ 
ods lead to consistent results. We now explain why we aim at 
using the R-L method instead of the O-P procedure. 

At this point of the analysis, it is necessa ry to briefly review 
the concept of the O-P method. We refer to iKriss et ^ (Il983l) 
for more details. The O-P method is a geometrical deprojection 
technique i nitially introduce d using the assumption of spherical 
symmetry (iKriss et al.ll983h . but it can be generalized to cl usters 
with axial symmetry (see, e.g.. lBuote & HumDhrevl2012ah . This 
method assumes that the cluster has an onion-like shell struc¬ 
ture, with a uniform emissivity in each of the concentric spheri¬ 
cal shells. The three-dimensional profile of the emissivity is then 
recovered from the surface brightness map in an iterative way 
from ring to ring, progressing from the outskirts to the center of 
the cluster. 

To carry out this iterative procedure, the O-P method requires 
the errors to propagate from ring to ring as the observed pro¬ 
jected profile is generated. When a quantity is overestimated in a 
radial bin, the same quantity needs to be underestimated in an ad¬ 
jacent bin, causing the me thod to produce fluct uations in the de¬ 
projected profile (see, e.g.. lAmeglio et al.l2007l) . Thus, the errors 
associated with each bin are correlated. The statistical uncertain¬ 
ties associated with the reconstructed three-dimensional profiles 
are usually estimated using Monte Carlo simulations (see, e.g., 
lBuotell2000ll . which by the same means smoothes the systematic 
errors due to the fluctuations in adjacent bins. 
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Fig. 3. Temperature profile obtained from tfie R-L reconstructed 
emissivity profile via Eq. ||6l witfi 7 = 1.19 + 0.04, compared to 
temperature profiles obtained using different techniques. Dark 
blue: profil e derived from combined X-ray an d thermal-SZ mea¬ 
surements dPlanck C ollabora tion et al.|l2013h using the ideal gas 
assumption (lEckert et al. 120131): light blue : spectral fitting of 
XMM observations ( Snowden et al.l J 008h: pink: s pectral fit¬ 
ting of Suzaku observations (iKawaharada et al.l2010h : red: mean 
value and errors on the reconstructed profile for y = 1.19; 
gray: mean value and errors on the reconstructed profile for 
y - 1.19 + 0.04. The mean and the error bars were obtained with 
a Monte Carlo method based on the R-L deprojection. The re¬ 
constructed temperature profile has been normalized to the data 
points from lEckert et al.l (1201 3h within a spherical shell delim¬ 
ited by the radii in the range [183; 2515] kpc. 


In the R-L method, the error propagation and the fluctua¬ 
tions due to the deprojection of the instrum ental noise are con¬ 
trolled by a regularization term (lLucvlll9^ . The choice of this 
regularization term is adapted to the data themselves and van¬ 
ishes as the deprojection converges. Eurthermore, simulations 
show that large-scale features are reproduc ed quickly, while 
small-scale structures are reco vered slowly dReblinskvl l2000t 
iPuchwein & BartelmanniUoObh . Therefore, the R-L method al¬ 
lows controlling th e fluctuation s that arise from deprojecting the 
instrumental noise (lLucvlll9^ . 

An important difference between the two methods is that the 
O-P method is nonlocal in the sense introduced above: to recon¬ 
struct the three-dimensional shell at a given radius, the results of 
the deprojection of all shells at larger radii are needed, and re¬ 
gions corresponding to different projected radii are mixed. The 
R-L deprojection method, however, is local in the sense that the 
two-dimensional profi le is independently deprojected along each 
line o f sight (see, e.g.. lReblinskvl20()0tlPuchwein & BartelmannI 
l2006h . such that the deprojected function is independently eval¬ 
uated at each image point and that this algorithm can be applied 
to data sets with incomplete coverage. However, the relations 
between different projected radii are controlled by the symme¬ 
try assumptions, and the information contained in different lines 
of sight may be mixed, even if the deprojection procedure is lo¬ 
cal. Eor instance, assuming spherical symmetry, the information 
collected at same projected radius is mixed for any line of sight. 

Nevertheless, as for the O-P method, it is possible to general¬ 
ize the R-L deproject ion method to intrinsically spheroiral clus¬ 
ter bo dies (see, e.g, iReblinsk^l2000t IPuchwein & BartelmannI 
I 2 OO 6 I) . In this formalism, the intrinsic shape of the projected 


body can be taken into account by a suitable (de)projection ker¬ 
nel in the R-L method. As the mixing between regions of dif¬ 
ferent projected radius is a consequence of the symmetries as¬ 
sumption and not of the deprojection procedure, testing different 
cluster symmetries should help us understand the effect of these 
mixings on the resulting profiles. 

This feature of local deprojection is most relevant for our 
study since the information provided by gravitational lensing 
measurements that we use here is also l ocal. In the SaWLens 
method developed bv iMerten et al.l (l2009l) . estimates of the lens¬ 
ing potential are obtained locally in all cells of a grid covering 
the region considered on the sky, and therefore, the observa¬ 
tions of weak gravitational lensing give information on the local 
galaxy density. 

3.3. Reconstructing the lensing potential 

The reconstructed projected potential ij/ is the projection along 
the line of sight of the gravitational potential rf). In our formal¬ 
ism, the gravitational pote ntial is expressed in terms of the recon¬ 
structed emissivity by (see lKonrad et al.l2013L for the derivation) 

jxirf , rj = . (7) 

The three-dimensional gravitational potential reconstructed this 
way is then simply projected along the line of sight. In this equa¬ 
tion, we use the bolometric emissivity {j^) and not the emissivity 
restricted to the energy band [0.4;2] keV. To recover the bolomet¬ 
ric emissivity from the energy-restricted emissivity (see Pig.fBi 
we used the t emperature profile (see Pig. ID and the thin-plasma 
code APEC (ISniith et alJl200lL fixing the metallicity to 0.3) to 
compute the ratio between the bolometric flux and the restricted 
flux for each radius. This factor was then used to recover the 
bolometric emissivity from the restricted emissivity. 

The lensing potential is recovered from observations of the 
distorted images of distant background galaxies. The distortion 
(shear and magnification) of these images is related to linear 
combinations of second derivatives of the lensing potential of 
a foreground cluster. This distortion is due to the gravitational 
tidal field of a deep potential well between the observer and the 
background galaxies. In the case of weak-lensing measurements, 
the distortion of individual galaxies is too weak to be observed, 
but the distortion averaged over sufficiently many neighboring 
galaxies is used to estimate the gravitational tidal field of the 
foreground galaxy cluster. In a large enough sample, the random 
intrinsic orientation of the background galaxies is expected to 
average out, and the ellipticity induced by lensing can be mea¬ 
sured. The lensing potential data used here were recovered with 
SaWLens in a nonparametric way from weak-lensing data. 

In Pig. m we show the reconstructed, two-dimensional po¬ 
tential assuming a polytropic index y — 1.19 ± 0.04. This profile 
is compared with the lensing potential obtained from weak lens¬ 
ing. At radii > 500 kpc our reconstructed potential profile agrees 
well with the reconstructed lensing profile within the error bars, 
even though the potential reconstructed from the X-ray emission 
appears to be slightly less curved. The discrepancy at small radii 
and the slight curvature change at larger radii may be caused by 
the lack of resolution in the weak-lensing measurement and by 
the other assumptions made (sphericity, polytropic and equilib¬ 
rium assumptions). These points are discussed in more detail in 
the next section. 

As we have mentioned in the introduction, studying the in¬ 
tracluster gas may help us understand the dark-matter distribu- 
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Fig. 4. Projected gra vitational potential re constructed using the 
method described in iKonrad et alJ (1201 3h compared to the po¬ 
tential recovered from weak gravitational lensing as a function 
of the projected radius s. The projected potential reconstructed 
from the X-ray emission has been obtained assuming a poly¬ 
tropic index y = 1.19 + 0.04 (see the text for details). The poten¬ 
tial recovered by lensing is shown in blue; the mean projected 
potential obtained from X-ray observations with y - 1.19 + 0.04 
is shown in red, the uncertainties of the latter are shown in gray. 
The mean and errors have been obtained using a Monte Carlo 
metho d to randomize the algorithm described in iKonrad et alJ 
(l2013h . The reconstructed lensing potential has been normalized 
to the lensing data within a circular shell region delimited by the 
radii in the range [92; 1734] kpc. 


tion in the cluster because the ICM fills the dark-matter poten¬ 
tial well. As a result of the locality of the R-L deprojection, our 
method opens a way for cluster potentials to complement the in¬ 
formation provided by lensing and joins two entirely different 
observables on the grounds of the gravitational potential, which 
provides a new tool for constraining the physical state of the gas 
locally. 


4. Discussion 

We used the R-L deprojection method to recover the three- 
dimensional gas distribution in Abell 1689 from its X-ray sur¬ 
face brightness. We then converted this gas distribution into a 
three-dimensional gravitational potential that we subsequently 
projected for comparison with the lensing potential. Using the 
R-L deprojection method is advantageous because it allows in¬ 
dependently deprojecting the X-ray observables along each line 
of sight as well as treating data sets with incomplete coverage. 
This therefore allows excluding regions where the equilibrium 
assumptions may not hold. As a result, a substantial amount of 
information becomes available by the direct comparison of the 
lensing potential with the projected potential reconstructed from 
X-ray emission (Fig. Ell. 

For example, the agreement between the two reconstructed 
potentials at radii from 500-2000 kpc indicates that the assump¬ 
tions underlying the derivation of the projected potential from 
X-ray observations, namely hydrostatic equilibrium, polytropic 
stratification, and the ideal-gas equation of state [Eqs. ([T]i-(Ell] 
are not substantially violated in the cluster Abell 1689 within this 
radial range. Nevertheless, we can observe a slight discrepancy 
between the potentials curvature: the potential recovered from 


lensing is flatter than the one reconstructed from the X-ray ob¬ 
servations. This flattening can be produced by the large resolu¬ 
tion limit present in the wea k-lensing measurements (~ 100 kpc, 
see, e.g., lBartelmannl2003L for a review), which can smooth the 
data so that the lensing potential recovered from weak-lensing 
measurements appears flatter than the potential recovered from 
X-ray observations. This discrepancy is expected to vanish if 
more potential observables are taken into account. For instance, 
a combined strong- and weak-lensing reconstruction is expected 
to yield a better match in the region close to the cluster center 
because the strong lensing is sensitive to this region. 


At small radii (R < 500 kpc), however, the two pro¬ 
jected potentials differ. Such a discrepancy in the cen¬ 
tral region of the cluster w a s also found by Peng et al.l 
(l2009f). [A ndersson & Mad eiskil (|20^. iMorandi et akl (1201 Ih . 


Miralda-Escude & Babull ( 1 99.51+ Lokas et ^ ( 2006h . and 


Serenoetal. (l2013l) . as mentioned in the introduction. This 


could be explained by a modification of the balance between 
the physical processes in the inner part of the cluster. It may 
indicate that some of the assumptions made should not be used 
in the central part of the cluster. 

Eor instance, the polytropic assumption with constant poly¬ 
tropic index may not be valid in the central part of cool-core 
clusters because the radiative cooling near the ce nter renders th e 
cooling time shorter than the Hubble time (e.g., lEabianl[T994l) . 
Cooling alters the polytropic equation by lowering the tempera¬ 
ture of the gas, which implies that the temperature profile does 
not increase toward the center as the gas density does, but in¬ 
stead decreases. Therefore, the assumption of a constant poly¬ 
tropic index is not expected to be correct at the cluster cen¬ 
tre. However, the assumption that the gas follows the polytropic 
equation seems to be valid in Abell 1689, with a fitted polytropic 
index of y= 1.19 ± 0.04 (Eig.|3]l. 


On the other hand, the assumption of hydrostatic equilibrium 
is likely valid in the central regions of cool-core clusters because 
of the deep gra vitational potential (see, for instance, Eig. 2 of 
iLau et al1l2009l) . This may n ot be the case in non-c ool-core clus¬ 
ters. As has been shown bv iMahdavi et al.l (1201 3h . for instance, 
departure from the equilibrium assumption is expected in un¬ 
relaxed clusters. The authors compared the lensing mass with 
the mass obtained under hydrostatic equilibrium assumptions for 
a sample of clusters at R < R 500 (~ 1.34 Mpc in the case of 
Abell 1689) and concluded that the departure from the equilib¬ 
rium assumption is different between relaxed and merging clus¬ 
ters: while for cool-core clusters, almost no bias is observed be¬ 
tween the lensing and the hydrostatic mass, a bias of 15-20% 
is present for non-cool-core clusters. Abell 1689 is known to 
experience a merging event aligned with the line of sight, and 
therefore, the gas is not expected to be in hydrostatic equilib¬ 
rium near the cluste r center (see also, e.g.. iNelson et al.l|2012l : 
iMorandi et al.ll20lTl) . This may contribute to the discrepancy at 
small radii observed in Eig. El 

Eurthermore, any departure from spherical symmetry could 
bias the estimate of the cluster mass, and if where the major axis 
of the cluster is preferentially aligned with the line of sight, it 
should lead to an overestimation of the lensing mass with respect 
to the X-ray mass. This is consistent with the conclusions of 
iLee & Sutol (120031 12 OO 4 I 1 . who showed that the dark matter pro¬ 
file in nonspherical clusters is always more sensitive to triaxiality 
than the gas distribution from X-ray measurements. This weak 
dependence of th e X-ray measurements on the cluster shape was 
also observed in iBuote & Humnhr^ (I2012alf5l . where the au- 
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thors computed the biases in the observables obtained by apply¬ 
ing a spherical O-P deprojection method on clusters with signif¬ 
icant triaxiality. They concluded that spherical averaging biases 
the observables by a small factor (< 1 %) with respect to a triax- 
ial deprojection. 

Therefore, the nonsphericity of the cluster could explain the 
discrepancy at small radii between the X-ray and lensing pro¬ 
jected potential profiles for Abell 1689, but interestingly, it does 
not seem to strongly affect the reconstruction of the potential in 
the region 500 - 2000 kpc. 

At large radii {R > R 500 ), simulations tend to show that the 
equilibrium assumptions are not valid any more due to mixing 
of the ICM wi th the infalling mater ial from the large-scale struc¬ 
ture (see, e.g. lReiprich et al.ll2013L for a review). The deviation 
from the hydrostatic equilibrium can be due to re sidual gas mo¬ 
tions or incomplete thermalization of the ICM (iMahdavi et alJ 
l2013h . Such phenomena are expected in the outskirts of clusters, 
where matter is accreted and substructures such as clurnps are 
more likely to be found (iNagai & Laull201ll : IVazza et alJl2013t 
IZhuravleva et al.ll2013h . so that our reconstruction method is not 
expected to be valid at radii larger than R 5 ooi~ 1.339 Mpc, for 
Abell 1689). However, as the lensing data used in our study 
do not extend to much larger radii than R 500 , we cannot test 
this statement. When we extrapolated the lensing data points to 
larger radii, we did not find significant differences between X- 
ray and lensing information in this radial range. 

Some systematic effects arise from the different assumptions 
we made in treating the X-ray observations. In Eq. lU, we as¬ 
sumed that the emissivity is dominated by bremsstrahlung emis¬ 
sion. This approximation represents a difference of at most 8% 
in the normalization compared to the normalization we would 
obtain taking into account the contribution of the emission lines 
for a metallicity of 0.3. Then, the conversion from 0.4-2 keV 
to the bolometric emissivity, computed assuming a temperature 
profile that follows Eq. 0 with y = 1.19, implies systematic 
errors of about 4% at most. The deviations from a flat field with 
the ROSAT satellite have been computed to intro duce system¬ 
atic eff ects on the order of 6% of the background in lEckert et al.l 
(l2012h . This value corresponds to the scatter in the residual, re¬ 
sulting from the fit of the surface brightness profiles of flat fields 
with a constant value. Einally, we estimated the effect of the lim¬ 
ited ROSAT PSE on the reconstructed profile by simulating a 
beta model with parameters consistent with those of Abell 1689 
and by convolving it with the PSE of ROSAT. While the con¬ 
volved profile is significantly flatter than the simulated profile 
inside ~300 kpc (by ~ 18%), beyond this radius the effect of the 
PSE is weak, and we estimated it to be at the level of ~ 3% of 
the observed signal. In general, these instrumental effects are not 
expected to change our conclusions qualitatively. Eurthermore, 
a combination of other cluster observables (such as SZ signal or 
galaxy kinematics) is expected to constrain the effects of these 
systematics and to help us build a consistent model for the pro¬ 
jected gravitational potential. 


5. Conclusion 

We applied a newly developed method (iKonrad et al.ll20lll for 
reconstructing the projected gravitational potential of the well- 
known galaxy cluster Abell 1689. We compared the X-ray pro¬ 
files obtained using the onion-peeling deprojection technique 


with the profiles from Richardson-Lucy deprojection. Eigures[T} 
|3]show the bremsstrahlung emissivity, electron density, and tem¬ 
perature profiles obtained using the two deprojection methods. 
All profiles are consistent within the error bars. 

Eigure |4] shows the projected gravitational potential recon- 
structed frorn X-ray data following the procedure described in 
iKonrad et al.l (1201 31) . This figure shows that the projected poten¬ 
tial obtained directly from X-ray observations agrees with the 
lensing potential obtained from weak-lensing measurements. 

Therefore, for sufficiently precise data sets, the R-L depro¬ 
jection method is expected to help us constrain the validity of 
the assumed equilibrium assumptions at all projected radii. We 
obtained these results using a spherical symmetry assumption. 
However, a possible bias coming from this assumption can also 
contribute to the discrepancy between the X-ray and lensing po¬ 
tentials at small radii (Eig. 4). To remove the degeneracy be¬ 
tween the different explanations (breakdown of hydrostatic equi¬ 
librium or projection effects), we are now working on generaliz¬ 
ing this deprojection method using a triaxial kernel. 

We provided a general proof of concept of potential reconstruc¬ 
tion combining X-ray data with weak-lensing measurements. 
The next step now consists of combining all observables pro¬ 
vided by galaxy clusters: strong and weak gravitational lensing. 
X-ray surface-brightness and temperature, galaxy kinematics, 
and the thermal Sunyaev-Zel’dovich effect into a joint construc¬ 
tion of one consistent model for the projected cluster potential. 
This combination will provide a substantial amount of informa¬ 
tion and will help us to quantify the validity of the assumptions 
made. If the assumptions are valid, the combination of all ob¬ 
servables will produce the best constrained galaxy-cluster po¬ 
tentials. 
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